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We discuss the implementation of the FONLL general-mass scheme for heavy quarks in deep- 
inelastic scattering in the FastKernel framework, used in the NNPDF series of global PDF analy- 
sis. We present the general features of FONLL and benchmark the accuracy of its implementation 
in FastKernel comparing with the Les Houches heavy quark benchmark tables. We then show 
preliminary results of the NNPDF2.1 analysis, in which heavy quark mass effects are included 
following the FONLL-A GM scheme. 
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1. FONLL in Mellin space: FastKernel 

The NNPDF series of PDF analysis [1-5] has used up to now the ZM-VFN scheme for the 
treatment of heavy flavours. In this contribution we review progress on the implementation of heavy 
quark mass effects in the NNPDF analysis and their implications on PDFs and LHC observables. 
We present preliminary results of the NNPDF2.1 set [7], which updates the global NNPDF2.0 fit 
with heavy quark mass effects as well as with the addition of the HI and ZEUS charm structure 
function data F2. 

Recently, the FONLL general-mass scheme, originally formulated for heavy quark hadropro- 
duction, was generalized to deep-inelastic scattering [6]. This scheme has several advantages over 
other existing GM schemes, such as S-ACOT [9], TR/MSTW08 [10] and BMSN [11]: it can be 
applied also in general hadronic processes, it can be formulated at any perturbative order, it allows 
the combination of different perturbative orders in the massless and massive computations, and it 
can be combined with various prescriptions for the treatment of subleading mass-suppressed terms 
near threshold, such as ;^-scaling or damping factors. In Ref. [6] explicit FONLL implementations 
are discussed at ff{a,) (FONLL-A) and at ff (a,^), either within an NLO (FONLL-B) or within a 
NNLO calculation (FONLL-C). 
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Figure 1: The F2 structure function for Q^ — 10 GeV^ computed at ff {as) accuracy (left plot) and i? (a^^) 
accuracy (right plot) in various GM schemes: FONLL, S-ACOT, MSTW08 and BMSN. The purely massless 
and massive results are also shown for comparison. 



The different versions of FONLL were compared to other commonly used GM schemes in the 
Les Houches heavy quark benchmark study [8]. Fig. |I] summarizes the result of this comparison. At 
(^(ttv), FONLL-A and S-ACOT are formally identical, while FONLL-B cannot be distinguished 
from the 0' (^aj) massive result at this value of Q^. In Fig. [l], the S-ACOT and FONLL-A curves 
have been obtained using a threshold damping factor prescription [6], the results with ;^-scaling 
are very close numerically. The MSTW08 result instead has been obtained using the Xi version 
of the ;t;— scaling prescription [8], which is rather different near threshold. The purely massless 
and massive results are also shown for illustration. It can be seen that all general-mass schemes 
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interpolate between the massless and massive calculations, and that there are significant differences 
between them, whose impact in PDF determinations is yet still to be systematically explored. 

The FONLL GM scheme has been now implemented in the NNPDF global analysis frame- 
work. As discussed in Refs. [2,5], for the PDF evolution and the computation of physical observ- 
ables the NNPDF analysis is based on the FastKernel framework, which allows a fast and accurate 
PDF evolution and evaluation of physical observables, including hadronic processes, without any 
K-factor approximation at all. The implementation of heavy quark effects in FastKernel has re- 
quired the computation of the MelUn transforms of the 6 {a^ massive DIS coefficient functions, 
both in the Neutral Current and in the Charged Current sector. These are not currently available in 
closed form; complete analytic expressions will be presented in Ref. [7], here we limit ourselves 
of presenting the results for the accuracy using the Les Houches heavy quark benchmark settings. 
We summarize the results of the comparison for F^ in Table |I|. In all cases the accuracy is well be- 
low the percent level, accurate enough for precision phenomenology. In Fig. g we compare the F^ 
structure function for various heavy quark schemes: ZM, FONLL-A-Damp and the FFN scheme 
as a function of ^ for different values of x. It can easily be seen how FONLL-A-Damp smoothly 
interpolates between the FFN scheme near threshold and the massless scheme at large Q^ . 
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Table 1: Results of the benchmark comparison for the F2c{x, Q^) structure function in the FONLL-A-Damp 
scheme for the FONLLdis code [6] and for the FastKernel framework. Results are provided at the benchmark 
kinematical points in x, Q^ [8]. 



2. Heavy quark mass effects in the NNPDF global analysis 



Here we present preliminary results based on a fit to the NNPDF2.0 dataset [5], supple- 
mented by charm structure function data F2 (x, 2^) from the HI and ZEUS experiments at HERA: 
NNPDF2.1. In Fig. ^ we perform a preliminary comparison of the singlet and gluon PDFs at the 
initial evolution scale Ql = 2 GeV^ in NNPDF2.1 and NNPDF2.0, normahzed to the NNPDF2.0 
central values. As expected, including heavy quark mass effects leads to a increase in the singlet 
at medium and small-x, as well as to a marked increase in the small-x gluon. Note however that 



Heavy quark mass effects in NNPDF 



J. Rojo 



U.1B 
0.16 
0.14 
0.12 
0.1 


FONLL-A-Damp — 
ZM-VFN 
FFN - - 


-■ 


^-"^y^ 




0.55 
0.5 

0.45 
0.4 

0.35 
0.3 


FONLL-A-Damp 

ZM-VFN 

FFN 


0.08 
0.06 
0.04 
0.02 



-^ 


'0 


, - ' 




0.25 
0.2 

0.15 
0.1 

0.05 



^-:><^-' 


1 


10 


100 


1 10 






Q^ [ GeV^ ] 




Q^ [ GeV^ ] 






x.io" 




x = 10"^ 


1.2 

1 

0.8 










2.5 

2 

1.5 




FONLL-A-Damp — 
ZM-VFN 
FFN - - 


-■ 


/ y^ ^ ' 




FONLL-A-Damp 

ZM-VFN 

FFN 


0.6 
0.4 




X 


,-' 




1 


y^y 


0.2 



^^ 


'<<-' 






0.5 



^^:>f^' 




10 


1 


DO 


1 10 






Q^ [ GeV^ ] 








0^ [ GeV^ ] 



Figure 2: The charm structure functions Fj (x, Q') as a function of Q^ for different values of x fromx = 10^^ 
to X = 10^^ in various heavy quark schemes, computed using the FastKernel method: FONLL-A-Damp, 
ZM-VFN and the FFN scheme. The PDFs and settings are identical to those of the Les Houches Heavy 
quark benchmark comparison [8]. 

differences are always within the PDF uncertainty bands. This implies that heavy quark mass ef- 
fects should modify the NNPDF2.0 predictions for LHC observables by 1-sigma or so at most. 
The implications of these results for LHC observables will be discussed in detail in Ref. [7]. This 
is to be compared with the CTEQ analysis [12], where for example a ~ 2.5-sigma variation on a^/ 
was obtained in the GM fit based on S-ACOT-;^ as compared to the ZM case. 

We also show in Fig. ^ the effects of heavy quark mass effects in the triplet and total valence 
PDFs: as expected their impact is completely negligible. Since processes which depend on valence 
PDFs, as well as those on the the medium and large-x gluon like Higgs production [13], are unaf- 
fected by heavy quark mass effects, the predictions for such processes obtained with NNPDF2.0 
will be very close to those of NNPDF2. L 
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Figure 3: Upper plots: The ratio of the preliminary NNPDF2. 1 Singlet (left plot) and gluon PDFs (right plot) 
to the respective NNPDF2.0 ones, at the initial evolution scale Qq—2 GeV^. Error bands are normaUzed to 
the NNPDF2.0 central value. Lower plots: same as above for the triplet Tj, and the total valence V. 
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